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We propose an experiment which can demonstrate quantum correlations in a physical scenario as discussed 
in the seminal work of Einstein, Podolsky and Rosen. Momentum-entangled massive particles are produced via 
the four-wave mixing process of two colliding Bose-Einstein condensates. The particles' quantum correlations 
can be shown in a double double-slit experiment or via ghost interference. 



I. Introduction 

Since the seminal works of Einstein, Podolsky and Rosen 
(EPR) [ 1 1 and Schrodinger |2| numerous experiments have 
demonstrated the counter-intuitive effects of quantum entan- 
glement, in particular the violation of local realism through 
Bell's inequality [3 1. Entanglement has been demonstrated for 
many physical systems such as photons 0^6], atoms 00, 
ions |9|[T0], and superconducting devices llTD . Variants of 
the EPR experiment have been realized, for example exploit- 
ing the analogy with quadrature phase operators |[T2l . But, 
until now, nobody was able to demonstrate the original EPR 
idea of an entangled state of freely moving massive particles 
in their external degrees of freedom, i.e. a state of the form 

J_oc dx I*>a I*+*o>b = J_ oo d /> ex P(i x 0P) \p)a |-P>b> where, x 
and p denote position and momentum, xq is a constant, and 
indices label the two particles. 

In this proposal, following the experimental approach of 
Ref. Ifl3l . we consider a Bose-Einstein condensate (BEC) of 
metastable helium-4 ( 4 He*). Via interactions with lasers, the 
particles are outcoupled from the trap, brought to collision, 
and fall on a detector. The collisions prepare atom pairs in a 
three-dimensional version of the EPR state with fixed abso- 
lute momenta. While the idealized EPR state would have to 
be created everywhere in space, our pairs are created in a fi- 
nite space volume. We propose to test the entanglement in a 
double double-slit experiment or via ghost interference. 

We start with a short review of the procedure described 
in Ref. Q3): 4 He* atoms of mass m - 6.646- 10 -27 kg are 
magnetically trapped and form a cigar-shaped BEC along the 
(horizontal) x direction. One employs two cr-polarized laser 
beams counterpropagating horizontally along the +x and —x 
directions and a 7r-polarized laser beam from the top along the 
-z direction fulfilling the Raman condition, all with a wave- 
length /1l — 1.083 |im. They bring the atoms into a mag- 
netically insensitive state and thereby induce a velocity kick 
in the horizontal +x directions and a kick upwards along +z- 
The recoil velocity of each kick is v rec = ^ 92mm/s, 
where h is Planck's constant. Thus, the three laser beams pro- 
duce a superposition of two counterpropagating matter waves 
of falling helium atoms, which subsequently scatter in a four- 
wave mixing process. In total, a fraction of about 5 % of all 



10 4 to 10 5 helium atoms collides and is scattered from the 
two condensates. While the relative velocity of two scattered 
atoms is 2 v lec , the velocity uncertainties can be obtained from 
the Gross-Pitaevskii equation. Assuming the trap parameters 
of Ref. [13], the BEC is elongated along +x and the uncertain- 
ties are anisotropic: Av x ^ 0.0044 v lec , Av v , ^ 0.091 v rec . 

The collisions are, to a very good approximation, of s- 
wave type, i.e. isotropic, and take place over a characteris- 
tic timescale of 150 |xs fl3l . Depending on the size of the 
condensate the collisions can produce momentum correlated 
particle pairs, lying on a shell in velocity space, whose origin 
is at v rec e z with radius v rec . Within quantum mechanical un- 
certainties, momentum conservation requires the two partners 
to find themselves being anticorrelated to each other in mo- 
mentum space. Most importantly, the isotropic nature of the 
s-wave scattering process gives rise to the superposition of all 
possible emission directions and thus to quantum mechanical 
entanglement in the external degrees of freedom of the two 
massive particles. 

The atomic de Broglie wavelength A dB associated with the 
recoil velocity is the same as the wavelength of the laser 
beams: A dB = » 1.083 urn. In Refs. |fT4l and 1731 a 

in Vrec 

model for the observed Hanbury-Brown Twiss (collinear) and 
back-to-back (BB) correlations of Ref. fPJl is developed. 



II. Double double-slit experiment 

Now consider a double double-slit arrangement as in Fig- 
ure[T] All particles hit by the lasers get a velocity kick of v rec in 
the ±x direction as well as upwards along +z. Let us consider 
only those particles which collided in such a way that they 
did not get any additional vertical velocity component and are 
moving along +y direction with velocity v rec after the collision 
(i.e. along +e y + e z with velocity "V2 v rec ). With gravity accel- 
eration g =i 9.81 m/s 2 and distance to the detector H = 0.5 m, 
the falling time is t - 328.8 ms. Having the coordinate origin 
O in the center of the initial condensate, the atoms pass the 
double-slit at a lateral position of y = L\ at some height -h 
and hit the detector at the maximally possible lateral distance 
y — L\ + L.2 = v rec t - 30.2 mm. 

We can ignore the effect of gravity by considering only the 
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FIG. 1 : Schematic of the double double-slit experiment (not drawn 
to scale). Pairs of atoms collide at height H and fall under gravity 
through the double slits onto the detector. See main text for details. 
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FIG. 2: Schematic of the double double-slit experiment (top view, 
not drawn to scale). Pairs of atoms, A and B, are emitted from source 
points r s within the BEC, pass double slits and arrive at detectors at 
positions r A and r B . See main text for details. 



top view of the experiment (Figure |2j. The source S of size 
S x x Sy x S z is emitting particle pairs A and B with veloc- 
ity Vrec- At a distance L\, there are two double-slits with slit 
separation d. At a further distance L^ an observation screen 
is located. There are four possible paths for getting a coinci- 
dence between atom detections on the left and right side at po- 
sitions Ta and Tb with x-coordinates x& and xb, respectively: 
{ai,bi}, {ai,b2}, {a2,bi}, and {a2,b2}. Here {a,,bt} means par- 
ticle A passed slit a, and particle B passed slit bk (z, k - 1,2). 

As discussed in Refs. [ 16- 18 ], there are two limiting cases 
in a double double-slit experiment: 

1. If the source size is very small, there is no momen- 
tum correlation between the particles. The momentum 
spread of each individual particle is large enough so that 
it can go through either slit at the same time and its 
partner does not carry enough which-path-information 
to identify through which of them it went. Thus, each 
particle forms a Young pattern independently of its part- 
ner. One effectively has two single-particle interference 
patterns unrelated to entanglement. The far field double 
double-slit (dds) two-particle pattern at the observation 
screens is a product of two independent one-particle 



patterns and is of the form 

lC S) ^A,XB)| 2 -COS 2 (^)cOS 2 ^ 

with df being the fringe distance. 



(1) 



2. If the source size is large (but still small compared to 
the slit distance), there is no one-particle interference 
pattern at either screen. The large source implies a 
small two-particle momentum uncertainty and therefore 
a high momentum correlation of the particle pairs. If 
one particle went through one slit, the other particle 
must have gone through the diagonally opposite one. 
However, since every particle is detected far behind a 
double-slit, the which-path information about its part- 
ner is erased. There is a superposition of two possibili- 
ties: particle A via upper slit & particle B via lower slit, 
and vice versa. The far field two-particle pattern cannot 
be factorized and has the form 



1, (dds), s,2 2/ X A-*B 

\lf/\ R '(X A ,X B )\ oc COS \7T 



(2) 



The Young fringes on, say, side B can be seen only con- 
ditionally on finding particles at a certain detector po- 
sition jca on the left side. Only by measuring coinci- 
dences, an interference pattern can be seen. Its maxi- 
mum on the right side is at the same x-position as the 
detector at the left side. 

Importantly, there is also a third regime, namely the one of 
very large source size. If the condensate is comparable to or 
larger than the slit distance, the two-particle interference pat- 
tern arising from the diagonal paths {ai,b2} & {a2,bi} is super- 
posed by a two-particle pattern originating from the horizontal 
paths {ai,bi) & {a2,b2). Therefore, detection of one particle 
at either slit does not imply any information about the slit the 
other particle goes through. This restores one-particle inter- 
ference and, due to complementarity, no genuine two-particle 
interference arises fi9l . 

To calculate the two-particle interference pattern on the ob- 
servation screen, we follow the treatment in Ref. [18|, where 
one integrates over point sources which emit two spherical 
waves without any (anti)correlation in momentum. Remark- 
ably, the anticorrelation and entanglement emerge naturally 
by integrating spherical waves of two particles emitted from 
the same position over a sufficiently large source area. This 
is in analogy to the Fourier transformation of a single parti- 
cle in phase space, giving rise to reduced momentum uncer- 
tainty as the source grows larger. Let us denote the possible 
path lengths from rs = (xs, ys, zs) to Ta = (xa, L\ + L2, 0) and 
r B = (xb, -Li -Li, 0) by L a . = ria^-alrA - and L bi = r s b,+b,r B 
with i = 1,2, as calculated by simple geometry. The (un- 
normalized) quantum mechanical amplitude for two entangled 
particles, emerging from rs, to land at points Ta and Tb is 



^SAB( r S> r A.rB) oc e 1 



-(£.,+**,) + e i 4 a "' +ib 2 ) 



(3) 
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It is a superposition of four equal-weight amplitudes cor- 
responding to the possible path combinations of particles 
A and B. We have omitted the one-over-distance dependence 
of the amplitude of the spherical waves because it is practi- 
cally constant in the far field and can be taken into the nor- 
malization factor. The phase of the amplitude of a collision to 
happen can be treated constant over the source volume if one 
neglects the expansion of the colliding BECs. Then, we can 
write down the quantum mechanical amplitude for two en- 
tangled particles, emerging from the whole source, to land at 
points Ta and Tb- It is a superposition of all possible emission 
points over the source volume V and is obtained via integra- 
tion over the whole condensate, possibly with some weighting 
function g(rs): 

^ s) (r A , r B ) oc I fff s dr s g(r 5 ) ifr SAB (r s , r A , r B ). (4) 

Summing up, we have the following conditions for a two- 
particle interference experiment: 

(I) The source must be sufficiently large to achieve well 
defined momentum correlation and wash out the single- 
particle interference pattern: 

^ « i (5) 
P L x 

The relative momentum spread — = must be 
small enough not to "illuminate" both slits. The source 
size S x implicitly influences Ap x . The larger the source 
size along x, the smaller Ap x becomes and the better the 
condition can be fulfilled. 

(II) The fringe distance df should be much larger than the 
detector resolution Sx: 

i,.*%>5*. (6) 
a 

Here we assume that one needs at least 5 pixels per os- 
cillation. 

(Ill) The source must be sufficiently small not to destroy the 
two-particle interference pattern: 

S x « d. (7) 

If the source size becomes comparable to or larger than 
the slit distance, two two-particle interference patterns 
wash each other out. 

Now we come to the parameter analysis. We take = 

I Vf sf - 0J Yf^ USD and S y = 10 Mm- This leaves s x as a 
free parameter, and we can write conditions (I), (II), and (III) 
in a single line: 

Li (i) _ (in) , (ID A dB L 2 

0.56io.mx — <k S x <§c d < . (8) 

d 5 6x 



S, = 25 urn, S, . = 1 nm S, = 50 Jim, S y = 10 |im 




FIG. 3: (Color online.) Two-particle probability distribution liA^'l 2 
for the double double-slit experiment for different source sizes. The 
slit distance d = 100 \im, source-slit distance h\ = 5 mm, and slit- 
detector distance L 2 = 25 mm, are kept fixed, resulting in a con- 
stant fringe distance df =s 271 \im. For a very small source (top 
left, S_ v = 25 (im) the momentum spread of each individual parti- 
cle is large and one obtains a product of two one-particle patterns 
of the form cos 2 (n^) cos 2 (n^). If the source is larger than the slit 
distance (bottom right, S x = 200 (i.m), two two-particle patterns to- 
gether wash out again into a factorizable pattern. An intermediate 
source size (top right, 5 , = 50|J.m) fulfills all conditions for two- 
particle interference and shows a distribution of the unfactorizable 
form cos 2 ^ ^' 8 ). 



An approximate solution for all our simultaneous constraints 
is d = 100 |J.m, L\ — 5 mm, L 2 = 25 mm. Here L\ + L2 ex- 
ploits the maximal possible lateral distance given by the dis- 
tance between BEC and detector of H — 0.5 m. The fringe 
distance becomes df ^271 \xm, which means that a fringe is 
resolved by only 4 to 5 pixels, given a detector resolution of 
Sx =t 60 [im. 

Figure [3]shows |(A^ S) | 2 for various source sizes S x . (The in- 
tegration was done over a two-dimensional rectangular source 
with size S y = 10 [im along y and constant weighting func- 
tion. Integration over y only marginally changes the pattern 
and so would an integration along 2.) While in the top left 
picture (S x = 25 |xm) condition (I) is violated, the bottom 
pictures (S x = 100 |xm and S x = 200 (i.m) violate condi- 
tion (III). The top right case (S x - 50 p.m) shows a con- 
ditional interference pattern with high genuine two-particle 
visibility. For this choice of parameters conditions ^ read 
28 \xm <K 50 \xm <K 100 |axn < 90 [im and are approxi- 
mately fulfilled. Decreasing or increasing the source size just 
a bit, lets us run into one or the other limitation. This shows 
that there is essentially no further freedom in any of the pa- 
rameters given typical experimental constraints (falling height 
H, detector resolution Sx). A double double-slit experiment 
therefore requires careful adjustment of the experimental pa- 
rameters. 
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FIG. 4: Schematic of the ghost interference experiment (top view, 
not drawn to scale). Pairs of atoms, A and B, are emitted from source 
points r s within the BEC and arrive at detectors at positions r A and 
Tb. Only atom A passes through a double slit. See main text for 
details. 

III. Ghost interference 

However, it turns out that we can circumvent condition (III), 
Ineq. Q, by removing one of the double-slits and using a 
ghost interference setup ETl . as shown in Figure |4] 

The possible path lengths from rs to r A = (xa,L\ + L2,Q) 
and Tb = (xg, ~L\ -L2,0) are abbreviated as L aj = rsa, + a,r A 
(z =1,2) and Ly, = FsFb. The quantum mechanical amplitude 
for two entangled particles in ghost interference (gh), emerg- 
ing from point r s , to land at points r A and r B is 

Or s , r A> r B ) - e 1 = ^ +Lb) + e 1 £ {L ^\ (9) 

The quantum mechanical amplitude for two entangled par- 
ticles, emerging from the whole source, to land at points 
r A and r B is again given by integration over all point sources: 

^(r A , r B ) oc I dr s g(r s ) ^(r s , r A , r B ). (10) 

In the ghost interference setup with only one double-slit, con- 
dition (III) is not necessary any longer. The remaining con- 
ditions (I) and (II) in (jS} can be easily fulfilled. Let us for 
instance choose d — 50 \im, L\ — 5 mm, L2 = 25 mm. 
Due to the smaller slit distance as compared to the double 
double-slit scenario, the fringe distance on side A becomes 
cfp = /Ids 7f - 542 \xm. The fringe distance on side B 
can be calculated via elementary geometrical considerations. 
It is d ( p = A dB ?ii±k ~ 758 lim _ Both can easily be re- 
solved with modern detectors. Conditions (I) and (II) read 
S x » 56 \xm, d — 50 \xm <K 83 \im, leaving the only con- 
straint S x » 56 |am. 

Figure [5] shows |iA A g'| 2 for different source sizes. One can 
see the transition from one -particle interference (S x = 50 [J.m, 
condition (I) not fulfilled) to almost perfect two-particle in- 
terference (S x = 400|J.m) of the form cos 2 (jT^jy s -). Note 
that condition (III), S x <sc d, is clearly violated for the larger 
source sizes. Moreover, in contrast to the double double-slit 
experiment, the maximum of the interference pattern on one 
side is opposite to the conditioning detection on the other side. 

Recently it was reported in Ref. [22 1 that matter waves from 
colliding BECs violate the Cauchy-Schwarz inequality, rul- 
ing out a description in terms of classical stochastic random 
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FIG. 5: (Color online.) Two-particle probability distribution I^b'I 2 
for the ghost interference experiment for different source sizes. The 
source-slit distance L t = 5 mm and slit-detector distance L 2 = 
25 mm are the same as in Figure [5] The slit distance d = 50 [im 
is smaller now, leading to larger fringe distances d^ A> = 542 \xm and 

d^ } = 758 p.m at side A and B, respectively. A very small source 
(top left, 5 x = 50 p.m) results in a product of two one-particle pat- 
terns. Large sources (bottom, S x = 200 [im and 400 p.m) show two- 
particle interference and produce a pattern of the unfactorizable form 
cosV^a). 



variable theories. The two-particle patterns discussed in the 
present work cannot be explained by classical correlations as 
soon as the two-particle visibility exceeds | l24ll25l . A strict 
proof that they must arise due to quantum entanglement can 
be made by employing a separability criterion using modular 
variables |23|. 



IV. Pair identification 

Finally, we come to an important issue, namely the problem 
of identifying coincidences. The pulse duration of the lasers is 
of the order of 500 ns, which is negligible. The traveling time 
of the recoiled atoms from one side of the condensate to the 
center, i.e. a distance of about 100 \im (more, if the condensate 
size is increased), is 1 ms. Therefore, all collisions certainly 
happen within At co h a 1 ms, most of them likely within a frac- 
tion of that time. Ref. [13] states that the time constant in the 
decay of the collision rate is 150 |xs. As discussed above, let 
us consider a pair of particles which collide into the +y and —y 
direction, with unchanged velocity component v rec along +z. 
Due to the velocity uncertainty we should assume that one 
particle has a velocity along z of v rec and its entangled partner 
has v rec + Av z . The times when they hit the detector plate are 
t Vicc ^ 328.8 ms and T Viec+ & Vz - 329.7 ms. This means that two 
particles which must experimentally be identified as a coin- 
cidence can have a time spread of AT pa ; r 0.9 ms, which is 
of the order of the collision time scale At co h. Two-particle in- 
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terference requires to measure coincidences of entangled part- 
ners, i.e. identification of the correct subensemble on side A 
conditional on detection on side B. Since an identification by 
arrival time seems to be impossible, it appears to be necessary 
to reduce the laser intensity such that per shot, on average, 
only a few pairs collide in a way that they can reach the max- 
imal lateral distance L\ + L2 on opposite sides of the detector. 
Taking into account the finite detection efficiency makes the 
pair identification experimentally extremely challenging. It 
may be advantageous to replace the double slit (in ghost inter- 
ference) by a grating to increase the possible count rates. 

V. Conclusion 

An experimental demonstration of the original EPR 
gedanken experiment using momentum entanglement be- 
tween pairs of two colliding BECs is within reach. While 
the problem of pair identification is very challenging, the con- 
straints of source size, time-of-flight distances, and detector 
resolution are manageable and let conditional two-particle in- 
terference in a double double-slit configuration or a ghost in- 
terference setup seem feasible. 
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